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We give a pedagogical introduction (suitable to upper level physics undergraduates) to the field 
of ultrarelativistic heavy ion collisions. We pay particular attention to our understanding of the 
thermodynamic and hydrodynamic properties of the matter created in heavy ion collisions at RHIC 
energies. 



PACS numbers: 

Finding the most fundamental blocks of matter has long been considered to be one of the fundamental problems 
within the physical sciences. 

Since a couple of centuries, we have established that the fundamental building blocks of matter as we know it are 
atoms. At about the same time, we realized atoms are in fact not fundamental, but are themselves made out of more 
basic blocks, protons, neutrons and electrons. 

Since the seventies, we know that protons and neutrons are made of constituents called quarks and gluons (or 
sometimes partons). The quest of better understanding the fundamental components of matter has then hit an 
obstacle. Actually, two obstacles. The first, obvious one, is that we found nothing new since then. Higher and 
higher colliding energies (by the uncertainity principle, the minimal size probed is roughly the inverse of the probing 
energy) have not yielded any observational result that would indicate that quarks are not pointlike but are themselves 
composite states. 

However, there is a different, more subtle snag: Since special relativity and quantum mechanics were combined, we 
know that, because of the equivalence of mass and energy (E = p 2 + mc 2 ) and the uncertainity principle {AEAt ~ 
ApAx r~j h), the number of particles is not fixed but subject to quantum fluctuations. 

This introduces a subtelty into the whole concept of "A is made of B" . Any system is, potentially, made up of 
an infinite number of particles, and the contributions of arbitrarily complicated multi-particle systems need to be 
included in any calculation (For a technical introduction to the topics described here, I recommend 

For electromagnetism (or, to be more precise, quantum Electrodynamics) this does not matter much: The quantum 
many body effects ("vacuum polarization" , shown in Fig. [Tj upper panel) reduce themselves to a slight decrease of the 
effective charge seen by charged particles with distance. In other words, when the two electric charges come closer 
together, their binding energy becomes slightly stronger. If they are really close, it could become strong enough that 
"two charges with a potential" is not anymore a good description of the system, so many particle-antiparticle pairs 
will appear from the vacuum that the series in number of particles will diverge. However, the distance for this to 
happen is so short (or the energy so large. This energy is called the Landau pole) that it simply lies beyond the 
regime where electrodynamics is a valid theory. 

Quarks and gluons, however, interact in a slightly different way, seen in the lower panel of the same figure. It all 
comes down to the fact that in electromagnetism there are only positive particles, negative particles, and neutral 
photons. Photons can create particle-anti particle pairs, but otherwise do not interact, and the fact that these pairs 
are massive means the fluctuations are automatically cut off (the vacuum needs to fluctuate by at least twice the 
electron mass to create such pairs, and this is unlikely) 

In QCD, however, quarks can come in three colors and antiquarks in three anti-colors (the way to make a neutral 
particle is either to mix a color with an anticolor, or three colors together). Gluons (the photons of QCD, also massless) 
carry a color and an anticolor. In other words, gluons are also charged and can interact with one another, and have 
more ways of interacting with one another than quarks, since there are 8 gluons (one for each color and different 
anti-color) and only 3 quarks. 
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FIG. 1: Inter-particle forces in electromagnetism and QCD 
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FIG. 3: Deep inelastic scattering 



Most importantly, the effect of these gluon loops on vacuum polarization is opposite to that of the quarks, anti- 
screening: For small distances charges between quarks are small, for large distances effective charges become larger 
(gluon vacuum fluctuestions make them increase), and for a distance of about lfm they diverge: The system becomes 
some not yet understood coherent collective system of many quarks and gluons. If the original color charges were 
infinitely heavy, one can intuitively see what happens in Fig. [2 Field lines interact and (usually) attract. So if the 
distance between them is big, they converge to form a tube. And therefore, by Gauss's law the force between two 
charges stops decreasing with distance, since the area of the tube remains constant. In fact, due to the anti-screening 
effect it increases. Since most quarks are not heavy, but actually ~ the electron mass, the tube can only be stretched 
until a pretty small (~ 1 fm, 10 -15 m) distance. Beyond it, the energy field will become strong enough to make new 
stable quarks out of the vacuum, and the string "breaks" into new color-neutral quark-antiquark pairs. The mass 
of protons is, to an acceptable approximation, nothing else but the uncertainity principle estimate from this size, 
h/(lfm* c),~ lGeV. Below these distances, strings are short enough that they "stretch out", and charge is small 
enough that perturbation theory applies. 

The above reasoning suggests that a good description of strongly interacting particles (called "hadrons" ) is that of 
"bags" of free quarks: As long as the quark stays inside the bag, it behaves as a free particle (no room for strings). 
But if it is kicked out of the bag, it develops a string, that ends up breaking into quark-antiquark pairs which end 
up forming color neutral objects, hadrons moving close together. As long as we do not worry about the details these 
hadrons (treat each correlated bunch of fast hadrons as a "jet" and assume its energy is that of the original quark) 
this picture allows us to calculate cross-sections for deep- inelastic scattering e — p collision experiments Q Fig [3] which 
have allowed us to test weakly coupled QCD thoroughly. The theory works to the extent we are sure it is true. We 
just can not see the fundamental degrees of freedom directly. 

So, would it be possible to produce "free quarks" that can travel and interact over a large volume? The previous 
discussion does leave a way out: If the distance between a bunch of protons and neutrons becomes comparable to the 
proton size, the bags "melt" and quarks behave as normal particles, interacting via a weak potential. 

One way to do this would be to compress "a lot of" nuclear matter to the point where the average separation 
between the nuclei would be less than the size of the proton. Quarks would then be free to jump from bag to bag, so 
the bags would, presumably simply disappear. Another way to see it is to realize that the such a compressed system 
becomes so hot that quarks will be able to overcome the string-like potential and travel from quark to quark. Such a 
picture is strongly suggestive of a phase transition, similar to the phase transition between ice and water (which occurs 
when molecules get enough average energy to escape the intermolecular forces which keep them in their position in 
the ice crystal). 

What is the phase transition temperature from the "hadron phase" to the "quark phase" ? A quantitative estimate 
for this is to compare the proton size (lfm) with the density of an ideal gas of 4 types of Fermions (up and down 
quarks, spin up and down) and 16 types of gluons(8 and 2 polarizations). The density is (In natural units where 
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FIG. 4: What we need to study thermodynamics and what we get in heavy ion collisions 
h = c = 1 and he = 0.197GeV fm) 

n = So (^« arfcs + ^ 9 ta 8 ) ~ 2.136T 3 (1) 

If for quark liberation the mean distance between quarks, n 1 / 3 ~ 1/m it leads to a temperature of ~ 152 MeV (a 
trillion or so degrees). More sophisticated models based on lattice QCD fj] put this temperature to be a surprisingly 
similar 150-200 MeV, leading to the suggestion that soon after the phase transition quarks become a nearly ideal gas 
describeable by perturbation theory techniques. A millionth of a second after the big bang, the whole universe was 
such a gas of quarks and gluons, called Quark GLuon Plasma (QGP). 

We do not, however, know much more about this state of matter. Is it really weakly coupled or do strong interactions 
persist? How does the transition from quark gluon matter to normal matter proceed? Is it really a phase transition, 
like ice and water, or is it a cross-over, like atoms and electromagnetic plasma? 

To explore these questions empirically is to create a chunk of quark-gluon matter in the laboratory. Fig. U shows 
the only way we can think of to do it: We collide two "large" nuclei at very high energy, and hope that some of this 
energy goes into heat. Enough heat to create a locally equilibrated liquid of quark matter, whose thermodynamic 
properties we then study. 

The trouble is that this is not the same as heating up matter in an "infinitely slowly heating oven" where thermo- 
dynamics applies automatically (Fig [4]). Even if the quark matter was created, it will soon expand, transform into 
hadrons, and decouple. Finding an unambiguous experimental sign that quark matter was in fact created is not a 
trivial matter. In fact, confirming experimentally quark matter was created is a still unresolved question. 

A plethora of experimental probes have been proposed to resolve it, ranging from quark chemistry (strangeness 
enhancement), to elecromagnetic signals (photons, dileptons), to the study of fluctuations (In context of susceptibilities 
diverging in phase transitions), to the behaviour of very high energyparticles ("jets") and heavy quark bound states. 
The interested reader can find a good introduction to these in |5j, [6( and the current experimental situation in 
0) H, H EH HH ■ In a sentence, we do not know enough and giving a balanced picture would take many more pages 
than allowed here. So I will devote the remainder of this work to explore an even a more basic question: 

Phase transitions and thermodynamic properties are relevant for systems in local thermal and chemical equilibrium. 
It is not at all clear that a short-lived, fast expanding system is in local equilibrium to the extent that thermodynamic 
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FIG. 5: Thermal particle fit to particle ratios at RHIC @, fill 03 



concepts such as phase transition and temperature will be relevant at all. 

In fact, the idea to describe strong interactions by statistical mechanics techniques, where the final abundance of 
particles obeys some kind of thermal distribution in terms of the particle's degeneracy g* and temperature (taking 
into account feed-down from resonances decays with branching ratios bj— »j) 
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This thermodynamic approach has a long and illustrious history [13, E3, E3, EH . It was using these techniques that 
[IH has managed to predict that hadronic gas had a phase transition even before the existance of quarks 1 Today, it is 
very well known that a fitted temperature, chemical potential and volume can describe the abundances of all strongly 
interacting particles produced in heavy ion collisions 0, EH (Fig- El for the simplest model, although extensions are 
possible, see (f|) The temperature fitted is very close to the phase transtion temperature at highest energies, and a 
universal freeze-out condition such as Energy/Baryon ~ 1 at lower energies. 

The interpretation of this finding is, however, hotly contested, with some people attributing it to more mundane 
processes (in multi-particle collisions phase space weights dominate over dynamics, and that gives you approximately 
statistical mechanics 13]) to very exotic ideas (strongly interacting particles are "born in equilibrium" [l7j . like 
Hawking radiation [18 1 from a black hole). Related to this discussion is the controversy of how small does system 
have to be before statistical models are valid. Qualitatitatively, even e + e~ — > hadrons looks statistical [13, EH, EH • 
Quantitatively, different groups give different conclusions as to how good the fit it [13, EH Here we evade this 
controversy by requiring that matter is not just "globally in equilibrium" (particles look like they have a temperature) 
but "locally in equilibrium" (particles have a temperature, pressure and density at each point. See Fig. [6] for the 
difference. This, invariably, leads us to consider not just temperature, but hydrodynamic flow. In fact, flow can be 
used as a tool to ascertain the thermodynamic and hydrodynamic properties of the system. 

Quantitatively, a locally equilibrated "fluid" can be distinguished from a far-from equilibrium "dust" of particles by 
using the Knudsen number [20(, K = l m f p /L where l m f p is the mean free path between two "microscopic" collisions 



1 The reasoningwent as follows [15l |: it was noted that the statistical degeneracy of strongly interacting particles increases with its mass, 
so not only these particles can have arbitrary mass, the more massive particles have greater statistical weight. It was reasoned that 
for heavy enough particles, distinguishing a "particle" from a "thermalized fireball" will be impossible, so we might as well consider 
all particles to be blobs of hot matter ("fireballs are made of fireballs"). This constrains the partitiion function to obey the recursion 
relation 
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It can be shown that this forces the density of states p r 
diverges 



to be exponential, p(m) ~ exp [/3h™] Hence, at T > ftV the partition function 
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FIG. 6: Global vs local equilibrium 



and L the "macroscopic" system size. It is clear that when K <C 1, each particle in the system interacts "many many 
times" with the other particles, and the system is, to a good approximation, a fluid where locally thermal equilibrium 
has been achieved. If, on the other hand, K > 1 the system is a "dust" , where particles to a good approximation 
ignore each other and statistical mechanics concepts such as temperature are inapplicable. 

As we know from [2(| , the mean free path for a weakly coupled gas is 1/ (an) where a the scattering cross-section 
between two particles and n the particle density, which for the massless gas limit is s — An where s is the entropy 

4 

density. In turn, in leading order perturbation theory a is related to the average momentum as a ~ , and by the 

massless gas formula (p) — 3T. Putting everything together, we get the following estimate for the Knudsen number 
0, [H and Eq. Q] 
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(3) 



remembering that for deconfinement T = 200 MeV, for a nucleus L ~ 10 fm, and for perturbation theory to be valid 
g <C 4-7T (say ~ 1 stretches it), we get a Knudsen number of 0.4, out of equilibrium by 40%. Hence, our hope of the 
system locally equilibrating are not really justified if the system is really perturbative. 

What if the system is not perturbative? Well, the very short answer is we dont know how to calculate things than, 
but we naively expect viscosity to be small since the system is hot and strongly coupled (lots of collisions, and "on 
average" each collision brings the system closer to equilibrium by Boltzmann's H-theorem [20| ) 

Could we do better than that? Well, if g > 1 we can not anymore use perturbation theory. From quantum 
mechanics, however, we can guess that the mean free path can not be smaller than the De Broglie wavelength (2~0 | 
(if it could, a thin layer of this material could function as a detector capable of breaking the uncertainly principle). 
Hence, provided (p) ~ 3T is still valid the Knudsen number would be something like 
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0.01 



(4) 



so the system would be very nearly thermalized (even for p — p collisions K ~ 0.3) Of course, the (p) ~ 3T assumption 
is completely unjustified in such a strongly coupled system, and the mean free path is not well defined since the 
number of particles will be extremely ill-defined. So this is at best a very rough estimate. 

Can we do any better? Recently, it was discovered that a theory "similar to QCD" (QCD with a large number of 
colors and 4 supersymmetries) is "dual" to a string theory defined on a 10 dimensional space with a string coupling 
inverse to the field theory coupling constant [22j (Note that such a theory is "conformal", so the coupling constant 
does not run. This conjecture, and it has not been proven yet, is known as the AdS/CFT correspondence). 

This means that an "infinitely strongly coupled theory" with infinite number of colors is equivalent to a weakly 
coupled string theory defined on a 5-dimensional space. The latter is simply classical general relativity and Hawking's 
formula for black hole entropy [l8l |. "Equivalent" simply means that it is possible to define a dictionary to convert 
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FIG. 7: Left panel: A dust (particles move independently of the initial shape) vs a fluid (particles move under forces determined 
by the initial density gradients). Right panel: The expected expansion of a fluid originating from an ellipsoidal hot region in a 
non-central collision. 



any Gauge-invariant observable in the quantum field theory into a corresponding observable in the string theory. The 
latter can be calculated since this theory is weakly coupled, and the answers will "magically" be the same. 

In particular, the entropy density becomes the "area" of a 5-dimensional black hole (calculable from Hawking 
formula) , and the viscosity becomes the speed at which small perturbations of the black hole dissipate (calculable in 
general relativity) . This led to the "famous" entropy and viscosity formula for an "infinitely strongly coupled theory" 
(of a very particular type) 23] f = 37- Plugged into our previous Knudsen number formulae this gives again a 
Knudsen number of 0.01 or so. So the system should certainly be a perfect fluid in the strongly coupled limit. 

Can all these ideas give experimentally measureable consequences? Well, we know what effect viscosity has on a 
fluid: It acts as a friction, transforming gradients of flow into heat [20| . Hence, a "perfect fluid" 's final momentum is 
determined only by initial gradients in density. In contrast, in a viscous fluid particles dont interact with each other, 
so gradients in density have no effect on the final kinematics. This can be see in Fig. [7] Viscous hydrodynamics 
should interpolate between these two extreme cases. 

The assumption of perfect thermal equilibrium at every instant is equivalent to assuming that in some frame the 
Energy- momentum tensor is diagonal, and pressure is isotropic (the same in all directions), so [201 ] 

/ e 

rpcomoving I J> 

^ ' p 

V p 

where the projector A^ = u^u u — g^ v projects a rank-2 tensor in the direction perpendicular to the flow u^. Here 
g^ u is the metric tensor, p the pressure and e the energy density. The latter are related by statistical mechanics 
(for the ideal massless gas discussed earlier p = e/3). This means that, for a perfect fluid, the conservation equation 
d„T A " y = becomes solvable from any set of initial conditions, since there are 5 equations (4 conservation equations 
and an equation of state) and 5 unknowns (u x ^ y ^ z ,p, e). These can be considered as the equations of motion of a 
collective system to "0-th order" in the Knudsen number, assuming the "macroscopic scale" is given by the flow 
gradient du. 

For viscous hydrodynamics, we have the additional equations to define shear and bulk viscosity ry, £ in terms of 
temperature and the energy momentum tensor is more complicated 

= T^ v \ equiUbrium + n MI/ (6) 



T£ neral = eu»u»+pA^ (5) 



= -77 ^A a ^d a u v + A av d a u^ - ^A^A^d'V) - Qd a u a A» v (!) 

but still solvable provided we can calculate 77, Q in terms of the temperature. Since 77, Q ~ l m fp (p) n [201 ] these can 
be considered the equations of motion to first order in the Knudsen number 2 (The complicated looking dissipative 



2 In fact, these equations have problems with relativstic causality, and have to be extended with terms to second order in the Knudsen 
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FIG. 8: Experimental results from V2 modelling of hydrodynamics, from [27| . N par t, the number of participants, measures the 
overlap area of the collision, see previous figure 

tensor 11^ is simply the traceless flow-perpendicular part of the gradient of the flow). Note that these equations 
are non-linear, so in general they can not be solved analytically but just integrated numerically from a set of initial 
conditions. 

The physical problem to which this formalism is to be applied is shown again in the right panel of Fig. [7] Initially, 
the system in a typical non-central collision, is, approximately, an ellipse, so its density transverse to the beam is 

e(r,0) = e o (r)(l + 2ecos2(^-fo)) (8) 

where 4>r is the reaction plane, seen in Fig. [JJ 

Integrating the equations shown earlier with these initial conditions 3 is involved and can not be done analytically. 
However, it can be readily seen that, if viscosity is small, the initial gradient in density will transform into a gradient 
in transverse flow vt 

u T {r,4>)=u T o{l + 2e v cos2((t>-(t>R)) (9) 

At the moment, we do not know how this fluid breaks up into particles. However, this break-up has to conserve 
momentum locally. This means that the average transverse momentum, as a function of <fi will also be something like 

( PT ) {<}>) = (pt) (l + 2v 2 cos 2(0 - <t> R )) (10) 

and this leads to a Fourier coefficient in the average pr binned by (4> — 4>r). This coefficient is known as v 2 [24]. v 2 is 
calculable within hydrodynamics, measurable by experiments, and highly sensitive to viscosity during the initial stage 
(Since hydrodynamic forces tend to make the fireball rounder, it is thought that v 2 stops forming when the system 
expands (24|. Note that the role of v 2 in the final stages of the fireball evolution is currently controversial). 

Comparing theoretical calculations to experimental measurments of v 2 0, H, B E3] j it can be seen that the system 
created at RHIC appears to be a very good fluid, well in the strongly coupled regime described earlier (Fig [5]) 
[2SHEH3]- I n f ac tj viscosity can not be too much bigger than r\ ~ s/4n, the lower bound value calculated through 
the AdS / CFT correspondence [23| . This has made a lot of people excited at the prospect that string theory can be 
used to describe the system created in heavy ion collisions. 



number. See [StJ for more details on this 
3 Initial conditions also need to give the density profile parallel the the beam. This is usually taken to be "Boost-invariant"', essentially 
assuming that the fluid flows apart in the beam direction with a Hubble type initial velocity v z = z/t. See 30] for a description 
justification of this initial condition. Note that even in the opposite ("Landau") limit Il4l where there is no longitudinal flow and nuclei 
"stick together" , a Bjorken-like flow develops very quickly, after a few fm. The applicability of these conditions, and the effect of them 
on the results, is currently not very well undestood 



8 



In this work we have concentrated on proving that a collective locally thermalized good fluid is indeed created in 
heavy ion collisions. What we did not say much about are the properties of this fluid: Is it really made of quarks 
and gluons? What is the nature of the transition between this fluid and normal hadronic matter? At what center of 
mass collision energy does such a "perfect fluid" form? How big does it really need to be before these collective effects 
arise? At what temperature, if any, does the system become weakly coupled and perturbation theory can be used? 

Other than lack of space, the main reason I did not mention these topics is that there is not much consensus about 
them. However, a lot of experimental activity is planned to try and clarify these questions. The focus is to vary as 
much as possible energysystem size, type of nucleus, collision geometry etc, to see if "something interesting" (eg a 
jump in V2, or multiplicity, or whatever) happens in some particular regime. 

The recently opened LHC accelerator plans to conduct very high energy heavy ion experiments, where we hope 
the temperature will be much bigger than the critical temperature. Further into the future, experiments with lower 
energy are planned to determine the critical threshold for production of quark matter. Two experiments of this type 
with a lot of Russian involvement are planned at FAIR in Germany [28| and NIC A in Dubna [29( . 

In conclusion, I heartily recommend Baikal school students to explore this subject further. It can be very theoretical 
yet includes heavy contact with experimental data. It draws insights from fields as far apart as thermodynamics and 
string theory. And, above all, it is very far from being thoroughly, or even approximately understood. In other words, 
we need your help! 

G.T. acknowledges the financial support received from the Helmholtz International Center for FAIR within the 
framework of the LOEWE program (Landesoffensive zur Entwicklung Wissenschaftlich-Okonomischcr Exzcllcnz) 
launched by the State of Hesse. GT thanks the organizing committee for the Baikal school for theoretical physics for 
inviting me to this great school, and all participants for questions, comments and discussions. 
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